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ABSTRACT 


The  nature  of  the  interaction  of  atmospheric  flow  with  a mountain 
range  of  finite  length  is  investigated.  An  inviscid,  adiabatic  primi- 
tive equation  model  is  used  to  simulate  steady,  vertically-unsheared 
flow  past  a mountain  range;  various  cases  of  mountain  height  and  width 
are  analyzed.  A similar  model  for  two-dimensional  Boussinesq  flow 
past  a mountain  on  an  f-plane  is  solved  analytically  following  Merkine 
(1975)  for  comparison  to  the  numerical  solution.  Results  indicate  that 
the  atmospheric  response  is  quite  similar  to  that  observed  in  the 
Boussinesq  model.  Ridging  over  the  mountains  and  responses  which  damp 
in  the  vertical  are  observed  in  both  models.  Differences  are  noted  in 
the  magnitude  of  the  response  in  the  vicinity  of  the  mountains  and  in 
the  position  of  the  downstream  trough.  In  addition,  the  effect  of 
finite  difference  model  grid  resolution  is  investigated  by  comparing 
simulations  based  on  differing  horizontal  and  vertical  resolution.  Re- 
sults indicate  that  for  larger  scale  mountains,  the  effect  of  poorer 
horizontal  resolution  is  significant;  vertical  resolution,  however,  is 
relatively  unimportant.  As  the  horizontal  scale  of  the  mountains  is 
reduced,  vertical  resolution  becomes  increasingly  more  important  so 
that  the  effects  of  both  horizontal  and  vertical  resolution  are  signi- 
ficant. 
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I.  INTRODUCTION 

For  many  years,  it  has  been  clear  that  mountain  barriers  have  impor- 
tant effects  on  cyclones  and  long  waves  as  well  as  local  effects  such 
as  mountain  waves.  The  proper  incorporation  of  topographic  features  and 
their  effects  in  numerical  models  is,  therefore,  a prerequisite  to 
accurate  forecasts. 

Most  numerical  weather  prediction  models  use  a form  of  topography 
which  is  linearly  smoothed  according  to  grid  size.  Linear  smoothing  is 
a very  serious  restriction  as  may  be  seen  in  the  hypothetical  case  posed 
by  Charney  (1967) . An  infinitesimally  narrow  mountain  ridge,  oriented 
along  a meridian  and  extending  to  the  top  of  the  atmosphere,  would  act 
as  an  absolute  barrier  to  the  flow.  Smoo tiling,  however,  reduces  the 
mountain  height,  and  therefore  its  effect,  to  zero. 

An  equally  significant  problem  is  that  the  dynamical  mechanisms  by 
which  mountains  influence  large  scale  flow  are  still  not  fully  under- 
stood. Earlier  investigations,  reviewed  by  Queney  (1973) , dealt  with 
linearized  flow  over  topography  and  were  primarily  concerned  with 
smaller  scale  features.  Recently,  Merkine  (1975)  considered  the  inter- 
action of  vertically-sheared,  inviscid,  steady  flow  with  an  infinitely 
long  ridge;  the  formulation  was  based  on  approximating  the  horizontal 
momentum  by  the  geostrophic  momentum.  Merkine  and  Kalnay-Rivas  (1976) 
applied  the  same  model  to  the  three-dimensional  flow  over  an  isolated 
topographic  feature. 

The  main  purpose  of  this  study  was  to  investigate  the  nature  of  the 
atmospheric  response  to  mountains  of  varying  widths  and  heights.  The 
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primitive  equation  model  developed  by  Monaco  and  Williams  (1975)  was 
used  with  a cyclic  lateral  boundary  condition  to  simulate  steady 
westerly  flow  past  a mountain  ridge  of  finite  length,  oriented  perpen- 
dicular to  the  mean  flow  in  the  middle  latitudes  of  the  northern 
hemisphere.  In  addition,  the  ability  of  the  finite  difference  model 
to  resolve  various  features  of  the  flow  is  considered.  A simplified 
case  with  features  quite  similar  to  those  of  the  model  atmosphere  was 
solved  analytically  in  a manner  similar  to  Merkine  (1975).  This  solu- 
tion was  then  compared  with  numerical  solutions  obtained  using  various 
horizontal  and  vertical  resolutions  in  the  model. 
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II . DESCRIPTION  OF  MODEL  AND  INITIALIZATION  PROCEDURES 
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The  two-level  general  circulation  model  used  in  this  research  has 
been  described  by  Monaco  and  Williams  (1975).  A brief  summary  of  the 
model  together  with  changes  made  during  this  research  is  presented  here. 


A.  MODEL  DESCRIPTION 

The  Monaco-Williams  model  is  a modified  version  of  the  UCLA  general 
circulation  model  described  in  detail  by  Arakawa  and  Mintz  (1974) . The 
atmosphere  is  treated  as  adiabatic  and  frictionless  with  no  sources  or 
si:iks.  The  equations  are  written  in  the  vertical  a-coordinate  system 
where: 


a = 


(2.1) 


7T  = 


- P 


T 


where  p is  the  pressure,  p^  is  the  constant  pressure  at  the  top  of 
the  model  atmosphere,  and  pg  is  the  surface  pressure.  The  vertical 
boundaries  of  the  model  are  the  earth's  surface  and  the  200  mb  constant 
pressure  level;  hence,  from  (2.1),  it  follows  that: 


a = 1 at  P = Ps  (2.2) 

a = 0 at  p = PT  = 200  mb 


• au 

The  boundary  condition,  0 = — = 0,  is  applied  at  the  upper  and  lower 
boundaries.  The  atmosphere  in  between  is  divided  into  layers  of  equal 


A 0.  The  structure  and  the  vertical  distribution  of  variables  is  shown 


in  figure  1 for  the  two-layer  version  of  the  model.  The  prognostic 


variables,  which  are  horizontal  velocity,  temperature,  and  surface  pre- 


sure,  are  governed  by  the  primitive  equations: 


a ,7T  a ,7tu  a .w  a .no  , rf  a 1 a 1,. 

3—  (—  u)  + 3jr( u)  + 3— ( U)  + 3—  ( U)  - [ + (V  rr U 3—  — ] ) ] 7TV 

3t  mn  35  n 9f|  m 3a  mn  mn  35  n 3r|  m 


+ — [37-  + 3r‘]  = F_  r 

n 35  dE,  mn  5 


(2.3a) 


a ,tt  ..a  ,ttu  , a ,itv  . a ir  a . irf  ..  a 1 a i., 

3t  mn  aC  n 3f|  m 3o  mn  mn  35  n 3r)  m 


TTra<f>  air,  it  „ 

+ —[3—  + oa  3—]  = — f , 
m 3n  3r|  mn  r| 


(2.3b) 


2—  (2L_  c T)  + c T)+  c t)  + (£— ) K -5—  (—  c 9)  = 

3t  mn  p 85  n p 813  m p p^  9o  mn  p 


,ir  , . u 3tt  v aiT,  , it  _ 

TOa  3-( — + — 33  + — -3-  + — Q , 

3t  mn  n 35  m 3n  mn 


(2.4) 


9 I11  \ . 3 U\  3 In  V\  3 _ n 

3t  mn  85(  n}  3n  m aa(mn)  ~ 0 ' 


(2.5) 


pa  = rt  , 


= — Tra 


Here,  — = a cos  <}>,—=  a ; a complete  listing  of  symbols  for  the  above 
m n 


equations  may  be  found  at  the  beginning  of  this  report.  The  model's 
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diagnostic  variables  include  the  vertical  a-velocity,  a,  and  geopoten- 


tial, <p  . The  vertical  a-velocity  is  given  by: 


7T 

mn 


a = 


i 

/ 


3 ,7TV.  , 

+ -5r(— ) I 
dn  m 


da 


3 .it  , 
° 3t  mn 


(2.8) 


Geopotentials  are  determined  in  a manner  which  maintains  the  integral 
property  of  the  vertically-integrated  pressure  gradient  and  conserves 
total  energy.  This  requires  combining  the  vertical  integral  of  the 
hydrostatic  equation  (2.9a)  with  an  alternate  form  of  the  hydrostatic 
equation  written  in  terms  of  potential  temperature  (2.9b) : 


1 

= y $ da  - tt  J' 


aada 


ii=-c9 

3pK  p 


(2.9a) 

(2.9b) 


The  horizontal  distribution  of  variables  (figure  2)  follows 
Arakawa's  "Scheme  C"  which  is  designed  to  simulate  accurately  the  geo- 
strophic  adjustment  process.  A detailed  discussion  of  this  scheme  and 
its  properties  is  given  by  Arakawa  and  Mintz  (1974) . In  this  research, 
the  latitudinal  spacing  of  the  TT  points  was  fixed  at  4°.  In  the 
longitudinal  direction,  two  different  grid  resolutions  were  considered: 
coarse  (AX  = 4.5°)  and  fine  (AX  = 1.5°) . Cyclic  continuity  with  wave 
number  8 was  assumed  in  the  ^-direction.  To  avoid  the  use  of  an  ex- 
tremely short  time  step  due  to  the  convergence  of  meridians  at  the 
poles,  the  ^-component  of  the  pressure  gradient  force  and  the  divergence 
in  the  ^-direction  are  zonally  smoothed  as  noted  by  Monaco  and  Williams 
(1975).  The  time  integration  proceeds  in  sequence  of  five  time  steps: 
one  Matsuno  forward/backward  step  followed  by  four  leapfrog  steps.  As 
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t 


I 
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TT  7T  • 

noted  by  Williams  (1976),  the  definitions  of  ( — ) and  ( — a)  at  u and  v 

mn  mn 


points  in  the  finite  difference  forms  of  the  zonal  and  meridional  momen- 
tum equations  were  changed  to  the  more  complex  forms  described  by 
Arakawa  and  Mintz  (1974) . 


B.  ANALYTIC  INITIALIZATION 

Initial  conditions  were  similar  to  those  described  by  Monaco  and 
Williams  (1975) . Haurwitz  (1940)  and  Neamtan  (1946)  obtained  the  stream 
function  solution  to  the  linearized  barotropic  vorticity  equation. 
Phillips  (1959)  used  this  solution  in  the  forcing  function  to  obtain 
the  geopotential  from  the  non-linear  balance  equation. 

Following  these  results,  an  atmosphere  in  solid  rotation  (no  vertical 
shear)  was  simulated  by  initializing  the  model  with  the  following 
equations : 


ip  = - B a sin  <p 


= 


i B a2  (2ft  + B) 


(2.10) 


u = 


1 9^  . 

-tt—  = Ba  cos  <p 

a dtp 


v = 


a cos  <p  9A 


where  <*>'  is  the  geopotential  perturbation  and  B is  a constant;  the 


other  terms  are  defined  at  the  beginning  of  this  report.  By  setting 


20  -1 

B = — sec  , a zonal  wind  distribution  which  varies  as  the  cosine  from 
a 


the  equator  (where  u = 20  m-sec  ^)  to  the  poles  (where  u = 0 m-sec  ^) 


was  initially  specified.  The  initial  temperature  distribution  was  accord- 
ing to  the  NACA  standard  atmosphere  as  defined  in  Haltiner  and  Martin 
(1957) . 
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C.  INTRODUCTION  OF  TOPOGRAPHY  INTO  THE  MODEL 


An  initial  problem  to  be  overcome  was  the  adjustment  of  the  model  to 
the  topography.  Adjustment  (or,  in  this  case,  steady  state)  was  assumed 
if  all  surface  pressure  tendencies  were  of  the  order  of  0.25  mb/hour  or 
less  and  fluctuations  in  certain  arbitrarily  selected  wind  components 
were  of  the  order  of  0.1  m-sec  Vhour  or  less.  A special  scheme  was  used 
to  introduce  the  topography  into  the  model  in  order  to  minimize  the  time 
required  to  reach  steady-state.  The  height  of  the  lower  boundary  was 
initially  set  to  zero  everywhere.  Heights  at  mountain  grid  points  were 
then  incremented  during  the  integration  as  a function  of  time  until  the 
desired  heights  were  reached.  The  equation  used  at  these  grid  points  is: 

M sin2  ) , t < t 
2tp  - p 

4 . (2.11) 

s 

M , t > t 

- P 

where  $s  is  the  height  of  the  mountain  at  time,  t;  M is  the  desired 
mountain  height;  and  t is  the  period  of  time  over  which  the  mountains 
are  "built".  A typical  time  series  of  the  surface  pressure  tendency  at 
.a  grid  point  in  the  vicinity  of  the  mountains  is  shown  in  figure  3 for 
an  integration  of  the  two-layer  fine  resolution  model  with  t = 16  hours. 
Note  that  the  largest  tendency  (~15  mb/hour)  is  observed  at  mid-period; 
this  corresponds  to  the  maximum  rate  of  change  of  the  mountain  height  de- 
fined by  (2.11).  With  one  exception,  it  appears  that  steady-state  has 
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III.  ANALYTIC  INVESTIGATION 


Merkine  (1975)  considered  a steady,  Boussinesq  flow  with  no  vertical 
shear  past  a large  mountain  ridge  in  a rotating,  stratified  atmosphere, 
bounded  above  by  a horizontal  rigid  lid.  He  was  able  to  analytically 
solve  the  potential  vorticity  equation  with  appropriate  boundary  condi- 
tions; he  was  also  able  to  derive  and  solve  an  analytic  equation  for  the 
v component  of  the  wind.  He  performed  this  analysis  on  an  f-plane  with 
infinite  boundaries  in  the  upstream  and  downstream  directions;  in  addi- 
tion, he  assumed  the  flow  to  be  balanced  in  the  vertical  and  cross- 
topography directions  (similar  to  Hoskins  and  Bretherton,  1972;  and 
Hoskins,  1974) . If  these  same  assumptions  are  made  together  with  a 
lateral  boundary  condition  of  cyclic  continuity,  an  analytic  solution 
can  be  derived  for  a simplified  case  which  is  quite  similar  to  the 
numerical  integrations  which  are  presented  in  this  report.  The  two- 
dimensional  potential  vorticity  equation  (independent  of  y)  for  inviscid, 
adiabatic,  Boussinesq  flow  on  an  f-plane  may  be  written: 


— [(— 
dt  u3x 


f)  2® 
’ 3z 


3v  39. 
3z  3x 


0 


(3.1) 


Consider  the  solution: 


3v  39 
3z  3x 


= constant 


(3.2) 


In  order  to  determine  the  constant,  equation  (3.2)  is  evaluated  away  from 
the  mountains  where  there  is  no  horizontal  or  vertical  shear  to  give: 

39 


3z 


= constant 


(3.3) 
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where  0s(z)  represents  the  mean  state  potential  temperature  at  level  z. 

1 9$ 

Using  the  geostrophic  relation,  v = — , and  the  hydrostatic  relation, 

ft  t dx 

o 9<f 

0 = 9 (z)  + — -r — , equation  (3.2)  may  be  written: 
s g dz 


1 a2*  xt2  ^ 1 a2*  s2$  ^ a2*  1 ,a2*  x2  „ 

f~2N  + “2  TT  TT  + “ 72  to5  = 0 

9x  f 9x  9z  9z  f 
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(3.4) 
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where  N = -g-  -s — . 

0 dz 
o 

In  order  to  non-dimensionalize  this  equation,  the  following  scale 
factors  (also  used  by  Merkine  (1975) ) are  introduced: 


$ = n2h2$' 


Z = Hz' 


NH  . 

X = — X 


The  non-dimensional  form  of  (3.4)  may  then  be  written: 


(3.5) 


92<J>'  92<j>'  92<j>'  92$' 


9x' 2 9x'2  9z'2 


9z 


1 2 


(3.6) 


or,  dropping  the  prime  notation: 


$ (1  + $ ) + $ - $ " 
XX  zz  zz  xz 


= 0 , 


(3.7) 


where  the  subscripts  x,  z indicate  differentiation  with  respect  to  a 


particular  variable  - e.g.,  0 = 


9$ 


x H * T^is  to  the  equation 

derived  by  Merkine  (1975) ; he  further  showed  that  this  equation  can  be 
transformed  into  Laplace's  equation  under  a non-linear  transformation. 
First,  equation  (3.7)  is  differentiated  with  respect  to  z , and  <J>  is 
replaced  by  ip  through  the  following  equation: 


tp  = $ + z - 1 . 

z 


(3.8) 
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Here,  \p  is  the  non-dimensional  stream  function  for  flow  in  the  x-z 


plane.  In  a steady-state,  the  non-dimensional  potential  temperature, 

$ + z , is  constant  along  a streamline.  The  new  equation  in  terms  of 

ip  can  be  written: 


ip  \p  2 + (\p  2 + 1)  ip  - 2\p  \p  \p  =0. 
xx  z x zz  x z xz 


Now  a transformation  of  variables  is  made  so  that  ip  is  the  new  verti- 
cal coordinate.  With  this  transformation  equation  (3.9)  becomes: 


z + z , , = 0 , 
xx  \p\p 


(3.10) 


where  z is  now  a dependent  variable  which  represents  the  height  of 
the  isentropic  surfaces.  It  is  convenient  to  define  the  new  variable 


Vi  = z-  ty+l. 


(3.11) 


In  terms  of  £ the  basic  equation  and  boundary  conditions  are: 


U + y. . = 0 

xx  ipty 


y = 0 at  ip  = 0 


(3.12) 


h _ ,1/2  , 
U = — f (S  x) 


where  f (S  x)  is  a function  specifying  the  height  of  the  lower 
boundary.  Since  the  numerical  model  invokes  cyclic  continuity  as  the 
lateral  boundary  condition,  let 


. , _ 1/ 2 , ,„l/2  . 

f (S  x)  = cos  (S  x) 


(3.13) 


_i/2  r 2Tnm  ...  . . , ...... 

where  S = — - — = — ; in  dimensional  form  (3.13)  becomes 


2ttx* 

cos  ( — £ — ) / where  the  (*)  denotes  dimensional  form.  Setting  y = 
1/2 

cos  (S  x)G(4>),  it  may  be  shown  that  (3.12)  reduces  to: 

G''Op)  - SG(tp)  = 0 


G (0)  = 0 


(3.14) 


G(-D=- 


which  has  as  a solution: 


G(tp)  = A sinh  (S1//2ip)  + B cosh  (S1//2iJ;) 


(3.15a) 


Equation  (3.15a)  is  evaluated  at  the  boundaries  to  determine  the  constants 


A and  B with  the  following  result:  A = 

the  solution  to  (3.12)  may  be  written: 


, B = 0 ; hence 


H sinh  S 


h sinh(S1/2ii>)  ,„l/2  . 

U = ~ 77  77o  cos(S  x)  * 

H sinh (S1/2) 


(3.15b) 


Making  use  of  the  transformation  rules  and  the  equations  of  motion, 
it  can  be  shown  that: 


(3.16) 


and  that  the  v component  is  determined  by  the  system: 


Vxx  + VW  = ° 
v = 0 at  ip  = 0 


(3.17) 


v,j,  = - £ S1/2  f*  (S1/2x)  at  ip  = 


2 


By  letting  f ' (S1,/2x)  = - sin(S1/,2x)  and  v = F(^)  sin(S1,/2x)  , this 
system  may  be  solved  in  the  same  manner  as  above;  hence. 


, h cosh  (S  ' ij/)  . ,„l/2  . 

v(x,<J0  ' 1/2  ~r"  3in(S  x)  • 

H sinh  S 


(3.18) 


Note  that 


vtoP  = v(x'0)  = - H " ;i/2 
H sinh  S 


sin(S  x)  , 


Vbottom  V^x'  U 1/2 

H tanh  S ' 


sin(S1//2x)  , 


are  the  solutions  at  the  upper  and  lower  boundaries.  The  fully  dimen- 
sional forms  of  v at  the  lower  and  upper  boundaries  may  be  written: 


V = 


Nh  ,2ttx*1 

— — f srn  (— ) , 

tanh(— ) 


(3.20) 


V = 


Nh  . ,2iTx* , 

<_> 

sinh ( fL  ) 


(3.21) 


The  asterisks  are  used  here  as  a reminder  that  x and  v are  in  their 

dimensional  forms;  v^*,  vt*  represent  values  at  the  lower  and  upper 

boundaries  respectively.  Considering  the  lower  boundary,  v * may  be 

b 

written  as  the  sum  of  an  infinite  series: 


*“*n  . ,2-rrx*. 

V - - L ~ 2TTNH.  Sin  W * 
n=i  tanh(f275T) 


(3.22) 


In  a finite  difference  model,  this  summation  must  be  truncated  after 


wave  number  nmax  = 2^  (the  maximum  wave  number  resolvable  by  the  grid) 


)r 


The  height  of  the  mountain  ridge,  h , is  specified  in  the  following 
manner.  In  the  dimensional  form,  zfa*  , the  height  of  the  lower 
boundary  may  be  written: 


zb*  = h cos  (- 


= h f(x*)  , 


where  f(x*)  is  a functional  representation  of  the  shape  of  the  topo- 
graphy varying  from  0 at  the  base  of  the  mountain  to  1 at  its  top. 

In  this  investigation,  only  mountains  symmetric  (in  the  east-west 
direction)  about  the  ridge  line  were  considered  (see  figure  4) . The 
exact  functional  form  of  the  mountain  may  be  written: 


= i-E 


- £TT  < X < £TT 


= 0 . 


eir  < x < it 


f (-^ — ) may  also  be  represented  by  the  following  Fourier-cosine  series: 


c ,nx*  'O  ,2nnx*, 

f (~g— ) = aq  + 2 ' An  c°s  < £ — ) , 


where: 


A = 2 ' 


e n=l 


. 2 (1-cos  (n£TT) ) 

A = r— r 

n 2 2 

n 77  e 


The  height  of  the  mountain  ridge  (at  x*  = 0)  may  then  be  written: 


h = h + h , 

o t—i  n 

n=l 


where : 


h --  h A 

o o 

h = h A 

n n 


hQ  does  not  contribute  to  the  v component;  therefore,  (3.26b)  is  used 


in  equation  (3.22)  to  derive  v * . 

b 
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In  order  to  analytically  approximate  the  conditions  present  in  the 


f'l 


numerical  model,  the  series  approximations  for  v * and  v * were 

t D 

truncated  after  wave  number  n = . For  example,  a 30-point  grid 

can  resolve  only  15  wave  numbers;  therefore,  the  series  approximations 
are  truncated  after  15  terms.  The  vfa*  analytic  profile  for  a 
9°  x 1.5  km  (width  x height)  mountain  computed  on  a 30-point  grid  span- 
ning one  cycle  of  wave  number  8 is  shown  in  figure  5 where  the  origin 
of  the  graph  is  placed  at  the  mountain  top.  Since  a 30-point  grid  on 
wave  number  8 results  in  1.5°  resolution,  the  mountain  is  carried  (by 
non-zero  z^*)  at  five  grid  points.  Superimposed  is  the  same  analytic 
case  computed  with  twice  the  resolution;  this  results  in  the  mountain 
being  represented  by  11  grid  points.  Note  that  the  two  profiles  at  the 
lower  boundary  are  nearly  coincident. 

In  figure  6 is  shown  the  v^*  analytic  profiles  for  a 3°  x 1.5  km 
mountain  computed  using  a 30-point  grid  versus  the  same  case  using  a 
60-point  grid.  In  the  former  case,  only  one  grid  point  represents  the 
mountain,  while  three  points  were  used  in  the  latter  case.  Qualitatively, 
the  two  profiles  exhibit  the  same  general  flow  characteristics;  however, 
the  fluctuations  observed  in  the  profile  of  the  30  grid  point  case 
differ  significantly  from  the  smooth  profile  of  the  60  grid  point  case. 
When  the  mountain  width  was  reduced  to  1.5°  on  the  60-point  grid  (there- 
by allowing  only  one  grid  point  to  represent  the  mountain) , similar 
fluctuations  were  observed.  These  results  seem  to  indicate  that  one  grid 
point  is  an  inadequate  representation  of  a mountain  with  respect  to  the 
quantitative  aspects  of  the  flow  at  the  lower  boundary. 

In  figure  7,  the  vfc*  analytic  profiles  for  a 3°  x 1.5  km  mountain 
computed  on  a 30-point  grid  versus  a 60-point  grid  are  shown.  The 
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similarity  between  the  two  profiles  indicates  that  the  poor  resolution 
has  little  effect  at  the  upper  boundary.  As  noted  above,  a grid  with 
1.5°  resolution  appears  to  be  sufficient  to  accurately  specify  the  flow 
response  to  a 9°  wide  mountain  at  both  the  upper  and  lower  boundaries. 

Therefore,  in  order  to  examine  the  characteristics  of  the  analytic  re- 
sponse, consider  figure  8 in  which  the  analytic  profiles  of  v^_*  and  v^* 
are  shown  for  a 9°  x 1.5  km  mountain  computed  on  a grid  with  1.5°  reso- 
lution. Both  profiles  are  antisymmetric  about  the  mountain  indicating 
ridging  centered  at  the  mountain  top.  A downstream  trough,  commonly 
associated  with  flow  over  mountains,  is  located  at  the  eastern  boundary. 

However,  in  this  case,  the  feature  is  purely  the  result  of  the  assumption 
of  cyclic  continuity.  Had  we  assumed  the  lateral  boundaries  to  be  in- 
finite, we  would  have  observed  downstream  a permanent  deflection  to  the 
south  (as  noted  by  Merkine  (1975))  because  of  the  assumption  of  a constant 
Coriolis  parameter.  Merkine  (1975)  also  noted  that  a variable  Coriolis 
force  would  generate  a trough  in  the  downstream  region  (results  from  the 
numerical  model  support  this  observation).  The  lower  boundary  response, 
as  expected,  is  greater  in  magnitude  than  the  upper  boundary  response. 

Note  that  the  lower  boundary  maxima,  however,  are  located  closer  to  the 
ridge.  This  would  seem  to  indicate  that  the  smaller  scale  response  is 
damped  more  rapidly  above  the  lower  boundary.  At  the  ridge  itself,  both 
vfa*  and  v * are  zero.  It  may  be  seen  from  (3.20)  and  (3.21)  that  varying 
the  mountain  height  affects  the  amplitude  of  the  response  only.  This 
can  be  seen  from  figure  9 in  which  the  v^*  profile  is  shown  for  a 9°  x 
1.5  km  mountain  vs.  a 9°  x 3.0  km  mountain.  Increasing  the  height  of  the 
mountain  by  a factor  of  two  causes  a uniform  increase  (by  a factor  of  two) 
in  the  magnitude  of  the  response  at  both  the  upper  and  lower  boundaries. 
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This  is  not  true  in  the  interior  where  non-linear  distortions  occur. 
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IV.  NUMERICAL  MODEL  RESULTS 

A number  of  integrations  were  performed  using  various  forms  of  the 
model.  In  each  of  the  experiments,  a mountain  range  of  uniform  height 
oriented  north-south  was  centered  in  the  grid.  The  model  was  initialized, 
as  earlier  described,  and  integrated  forward  in  time  until  a steady-state 
was  achieved.  In  the  following  discussion  again,  fine  resolution  refers 
to  1.5°  longitudinal  grid  spacing  and  coarse  resolution  refers  to  4.5° 
grid  spacing. 

A.  BROAD  MOUNTAIN  RANGE  CASE 

A mid-latitude  mountain  range,  9°  x 1.5  km  (width  x height),  extend- 
ing from  34°N  to  66°N  was  the  first  case  considered.  Results  from  the 
various  integrations  are  as  follows: 

I 

1.  Six  Layer  Fine  Resolution  Model  Results 

The  wind  field  from  level  six  (a  = .91667)  is  shown  in  figure  10. 

Significant  features  include  ridging  over  the  mountains  and  weak  trough- 

ing  downstream.  Upstream,  the  flow  is  deflected  northward  by  the 

mountains;  downstream,  the  deflection  is  southward,  becoming  weakly  north- 

* 

ward  east  of  the  trough.  At  a distance  from  the  mountains,  both  up- 
stream and  downstream,  little  deflection  is  apparent.  Deflection  of  the 
flow  is  also  observed  one  grid  point  (~445  km)  to  the  north  and  south  of 
the  mountains;  this  distance  is  of  the  order  of  the  Rossby  radius  of 
deformation.  In  figure  11  is  shown  the  wind  field  at  level  one  (a  = 

.08333).  In  general,  the  effect  of  the  mountains  on  flow  at  this  level 
appears  to  be  minimal.  Consider  the  v profile  for  level  six  (a  = .91667) 
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at  44°N  shown  in  figure  12.  This  latitude  was  chosen  because  it  crosses 


the  mountain  near  the  center  and,  more  importantly,  because  it  was  re- 
presentative of  the  response  at  the  other  latitudes  which  cross  the 
mountains.  Superimposed  on  the  v-profile  is  the  analytic  solution,  v^*, 
at  a = 1.0  for  the  similar  analytic  case.  There  is  significant  similar- 
ity between  the  two  curves.  Both  profiles  indicate  ridging  at  the 
mountain  top  with  the  only  apparent  upwind  difference  being  one  of  magni- 
tude (the  model  showing  a greater  northward  deflection) . Downstream, 
as  in  the  analytic  case,  the  southward  maximum  is  in  the  same  position, 
relative  to  the  mountain  top,  as  the  northward  maximum;  however,  its 
magnitude  is  about  twice  that  of  the  northward  maximum.  Also,  downstream, 
the  model  profile  indicates  the  trough  position  to  be  about  75°  (longi- 
tude) east  of  the  mountains;  in  contrast,  the  analytic  solution  shows 
the  trough  at  the  channel's  eastern  boundary.  This  difference  is  to  be 
expected  because  the  Boussinesq  approximation  and  the  f-plane  and  rigid 
lid  assumptions  (present  in  the  analytic  case)  are  not  made  in  the 
numerical  model.  In  figure  13,  the  v profiles  for  levels  two,  four  and 
six  at  44°N  are  shown.  Noting  the  positions  of  the  northward  and  south- 
ward maxima  relative  to  height,  a systematic  westward  tilt  in  the  verti- 
cal is  observed;  this  tilt  is  not  observed  in  the  analytic  solution. 

It  is  felt  that  this  tilt  could  be  the  result  of  the  S-effect  and  other 
effects  not  included  in  the  analytic  solution.  As  noted  in  the  analytic 
case,  the  response  should  damp  with  height — this  is  especially  apparent 
in  the  magnitude  of  the  southward  maximum.  In  figure  14  is  shown  the 
analytic  solution  at  0 = 0 versus  the  model  profiles  for  a = .08333  and 
0 = .25.  Note  that  in  both  model  profiles  the  magnitude  of  the  response 
is  similar  to  that  of  the  analytic  case,  but  the  horizontal  scale  is 
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smaller.  In  addition,  the  systematic  westward  tilt  with  height  is  ob- 
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served  in  the  two  model  profiles. 

2.  Two  Layer  Fine  Resolution  Model  Results 

The  same  case  was  computed  using  the  two  layer  model  to  examine 
the  effect  of  vertical  resolution  on  the  numerical  simulation  of  the 
flow.  The  a = .75  velocity  field  was  compared  with  its  counterpart  from 
the  six  layer  model  and  only  minor  differences  (in  the  vicinity  of  the 
mountains)  were  observed.  In  figure  15,  the  v-profiles  for  a = .75  at 
44°N  are  shown  for  these  two  models.  There  is  little  difference  between 
the  two  curves,  both  showing  a ridge  centered  one  grid  point  west  of  the 
mountains  and  both  showing  weak  troughing  in  the  same  general  area 
downstream.  Some  differences  exist  in  the  vicinity  of  the  mountains, 
but,  in  general,  the  two  layer  model's  response  appears  to  be  quite 
similar  to  that  of  the  six  layer  model.  In  figure  16,  the  v-profiles  at 
a = .75  and  0 = .25  from  the  two  layer  model  are  superimposed.  As  in 
the  case  of  the  six  layer  model,  a westward  tilt  is  again  observed. 

3.  Coarse  Resolution  Model  Results 

In  order  to  investigate  the  effect  of  horizontal  grid  resolution, 
the  above  integrations  were  repeated  using  the  six-layer  and  two-layer 
coarse  resolution  models.  The  velocity  fields  at  0 = .25  and  0 = .75 
produced  by  these  two  models  were  quite  similar.  In  figure  17  is  shown 
the  six-layer  model  field  at  O = .75.  In  general,  it  shows  a deflection 
of  the  flow  similar  to  the  six-layer  fine  resolution  model,  but  differ- 
ences were  noted  particularly  in  the  vicinity  of  the  mountains.  In 
figure  18  are  shown  the  v-profiles  for  the  six-layer  fine  resolution  model 
versus  the  six-layer  coarse  resolution  model.  Note  that  the  magnitudes 
are  about  the  same  but  that  the  scale  of  the  response  is  less  in  the 
fine  resolution  model.  This  would  appear  to  emphasize  the  importance  of 
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horizontal  grid  resolution  in  the  simulation  of  finer  scale  features 


i 
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of  the  flow  over  topography. 

Consider  figure  19  in  which  the  lower  level  (0  = .75)  v-profiles 
for  the  two-layer  coarse  resolution  model  are  shown  for  three  different 
mountain  height  cases  (1.5,  2.4,  and  3.0  km).  The  shapes  of  these 
three  profiles  are  quite  similar  with  magnitudes  increasing  with  in- 
creasing mountain  height.  Clearly,  the  response  is  directly  proportional 
to  mountain  height  as  it  is  in  the  analytic  case  at  the  boundaries.  In 
addition,  although  the  response  at  the  mountain  top  is  non-zero,  it  is 
approximately  independent  of  mountain  height.  Hence,  at  least  some  of 
the  gross  aspects  of  the  analytic  solution  are  apparent  in  the  coarse 
resolution  simulations. 

B.  NARROW  MOUNTAIN  RANGE  CASE 

Next,  a narrower  mountain  range  (3C  wide)  extending  from  34°N  to 
66°N  was  considered;  mountain  height  was  again  set  at  1.5  km.  Since 
the  distance  between  grid  points  in  the  ^-direction  is  1.5°  in  the  fine 
resolution  model,  this  allows  only  one  grid  point  at  each  latitude  to 
represent  the  mountains. 

1.  Six  Layer  Fine  Resolution  Model  Results 

The  wind  field  from  level  six  (0  = .91667)  is  shown  in  figure 
20.  As  in  the  wider  mountain  case,  there  is  ridging  over  the  mountains 
and  troughing  downstream.  However,  the  magnitude  of  the  deflection  in 
the  vicinity  of  the  mountains  appears  to  be  much  greater  in  this  case. 
Note  that  one  grid  point  (1.5°  longitude)  west  of  the  mountains,  the 
flow  is  almost  directly  from  the  south  the  whole  length  of  the  mountain 
range.  Significant  deflection  of  the  flow  is  observed  one  grid  point 
(~445  km)  north  and  south  of  the  mountains  as  was  true  in  the  wide 
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mountain  case.  The  effect  of  the  mountains  on  the  upper  level  flow  at 
a = 0.08333  appeared  to  be  minimal  (as  in  the  wide  mountain  case).  In 
figure  21  is  shown  the  v-profile  for  level  six  (a  = 0.91667)  at  44°N 
superimposed  on  the  analytic  solution  at  a = 1.0  for  the  same  case.  The 
strong  northward  deflection  observed  in  figure  20  is  seen  here  at  grid 
point  14  as  a strong  northward  maximum  exceeding,  in  magnitude,  the 
analytic  solution  by  a factor  of  four.  While  the  downstream  southward 
maximum  is  not  coincident  with  its  analytic  counterpart,  it  has  the  same 
magnitude  as  the  northward  maximum.  This  symmetric  behavior  is  predicted 
by  the  analytic  solution  and  it  differs  from  the  result  in  the  wide 
mountain  case  where  the  magnitude  of  the  southward  maximum  was  about 
twice  that  of  the  northward  maximum. 

The  position  of  the  downstream  trough  is  closer  (by  3°  longitude) 
to  the  mountains  in  this  case.  Note  the  similarity  between  the  fluctua- 
tions observed  in  the  analytic  solution  and  those  appearing  in  the  up- 
stream v-profile.  It  would  appear  that  these  oscillations  are  a result 
of  allowing  only  one  grid  point  to  represent  the  mountains  interacting 
with  a mean  flow  from  the  west.  These  oscillations  are  not  present  in 
the  model's  downstream  profile.  In  figure  12  are  shown  the  v-profiles 
for  levels  two,  four  and  six.  Noting  the  position  of  the  southward  maxi- 
mum with  respect  to  height,  there  again  appears  to  be  a systematic  west- 
ward tilt  with  height;  however,  it  is  not  as  well  defined  here  as  in 
the  wider  mountain  case  (figure  15) . In  figure  23  is  shown  the  analytic 
solution  at  a = .08333.  As  in  the  wide  mountain  case,  the  amplitudes  of 
the  two  responses  are  similar,  but  the  horizontal  scale  of  the  model  re- 
sponse is  much  smaller. 
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2.  Two  Layer  Fine  Resolution  Model  Results 


An  integration  for  the  same  case  was  performed  using  this  ver- 
sion of  the  model.  Comparing  the  resulting  wind  field  at  0 = .75  with 
the  six  layer  model  results,  significant  differences  were  observed  in 
the  vicinity  of  the  mountains  (see  figure  24) . In  the  six- layer  model 
field,  a well-defined,  narrow  trough  broadening  toward  the  south  is  ob- 
served immediately  east  of  the  mountains;  however,  there  is  no  trough 
in  this  area  in  the  wind  field  of  the  two  layer  model.  Differences 
were  noted  in  the  magnitudes  of  the  wind  velocity  downstream.  Fur- 
ther upstream  and  downstream,  the  two  models  give  similar  results  at 
this  level.  A comparison  of  the  fields  at  0 = .25  showed  little  differ- 
ence between  the  two  models.  In  figure  25  are  superimposed  the  v 
profiles  from  both  models  for  a = .75  at  44°N.  As  can  be  seen,  little 
differences  exist  away  from  the  mountains;  however,  immediately  upstream, 
the  oscillations  in  the  six-layer  model  profile  appear  to  be  of  greater 
amplitude,  and,  downstream,  the  six- layer  model  profile  shows  the 
presence  of  a strong  trough  which  is  absent  in  the  two-layer  model  pro- 
file. In  figure  26,  the  v-profiles  from  both  models  at  0 = .25  (also 
at  44°N)  are  shown;  the  two  are  similar  to  the  extent  of  having  a prin- 
cipal trough  and  ridge,  but  sizable  phase  differences  are  seen.  These 
results  indicate  that  the  vertical  resolution  is  important  in  this  case. 

During  the  latter  stages  of  this  study,  a relatively  minor  coding 

error  was  discovered  in  the  computer  program;  this  error  had  been 

present  throughout  the  period  of  this  research.  As  described  by  Monaco 

and  Williams  (1975) , a pressure-area-weighted  form  of  the  horizontal 
TTU  7TV 

wind  components  ( — , — ) is  forecast  by  the  model;  u and  v are  computed 
mn  mn 

TT 

after  each  time  step  by  dividing  out  — . The  error  resulted  in  a slight 


overestimation  of  this  term  in  that  section  of  the  coding  responsible 
for  removing  this  weighting  from  the  v component  after  each  leapfrog 


time  step.  A test  case  was  recomputed  using  the  corrected  version  of 
the  two-layer  coarse  resolution  model.  Almost  no  change  was  observed 
in  the  wind  fields  for  the  o = .25  and  a = .75  levels.  Although  slight 
changes  in  the  v-profiles  were  observed,  the  overall  pattern  remained 
the  same.  It  was,  therefore,  concluded  that  the  error  in  coding  did 
not  significantly  affect  the  nature  of  the  response  depicted  in  these 
results. 


V.  CONCLUSIONS 
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From  the  results  of  the  numerical  simulations,  it  is  evident  that, 
in  general,  the  mountains  interact  with  the  atmospheric  flow  in  such  a 
way  as  to  induce  ridging  over  the  mountains  and  troughing  on  the  down- 
stream side  at  lower  levels.  With  the  exception  of  the  troughing  the 
numerical  solution  appeared,  in  general,  to  agree  with  the  analytic 
result  from  the  Boussinesq  model.  The  presence  of  the  downstream  trough 
was  due  to  the  fact  that  the  numerical  model  has  a variable  Coriolis 
parameter  and  a non-rigid  upper  boundary;  in  contrast,  the  analytic 
model  considered  Boussinesq  flow  on  an  f-plane  bounded  above  by  a non- 
rigid  lid. 

Several  other  similarities  were  observed  between  the  analytic 
response  and  the  numerical  response.  As  in  the  analytic  case,  smaller 
scale  features  appeared  to  damp  rapidly  with  height.  Varying  the 
mountain  height  was, in  general,  observed  to  affect  the  amplitude  of  the 
numerical  response  in  the  same  manner  as  in  the  analytic  case.  Reduction 
of  the  width  of  the  mountains  reduced  the  horizontal  scale  of  the 
numerical  response,  but  in  contrast  to  the  analytic  solution,  it  also 
increased  the  magnitudes  of  the  northward  and  southward  maxima.  It  was 
initially  felt  that  the  strong  northward  maximum  along  the  western  side 
of  the  3°  x 1.5  km  mountains  might  be  related  to  air  flowing  around 
rather  than  over  the  mountains.  A detailed  examination,  however, 
seemed  to  indicate  that  this  was  not  the  case. 

The  numerical  solution,  in  every  case,  displayed  a systematic  west- 
ward tilt  with  height  which  was  not  observed  in  the  analytic  solution. 
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The  precise  reason  for  this  tilt  was  not  determined;  however,  it  appears 
to  be  the  result  of  the  3-effect,  compressibility  and  non-rigid  upper 
boundary  in  the  numerical  model. 

The  effect  of  vertical  resolution  appeared  to  be  minimal  in  the 
case  of  the  broader  mountains;  but  significant  effects  were  observed 
when  the  mountain  width  was  reduced.  With  respect  to  horizontal  resolu- 
tion, the  coarse  resolution  model  appeared  to  be  able  to  reproduce  some 
of  the  larger  scale  features  observed  in  fine  resolution  solutions; 
however,  differences  were  noted  between  the  two  models  with  respect  to 
finer  scale  features  of  the  flow. 

Early  in  this  research,  an  unsuccessful  attempt  was  made  to  simulate 
conditions  in  which  air  flows  around  rather  than  over  the  mountains. 

From  these  results,  and  those  of  Merkine  (1975),  it  appears  necessary 
to  further  reduce  the  horizontal  scale  of  the  mountains.  These  results 
also  indicate  that  finer  horizontal  and  vertical  resolution  would  be 
necessary  to  accurately  simulate  the  atmospheric  response  to  this 
narrower  mountain  range.  With  the  high  price  paid  in  computing  time 
for  increased  resolution,  this  appears  to  be  a definite  argument  for 
the  use  of  a nested  grid  in  treating  atmospheric  response  to  mountains. 

No  problems  were  encountered  in  treating  unsmoothed  topography  m 
the  numerical  model.  In  addition,  the  scheme  used  to  introduce  the 
mountains  into  the  model  behaved  quite  well.  Additional  attention 
should  be  given  to  the  possibility  of  reducing  the  period  of  time  neces- 
sary to  incorporate  the  mountains.  If  this  can  be  done,  a similar 
scheme  might  be  used  to  introduce  smaller  scale  topographical  features 
into  operational  numerical  models. 


Finally,  this  investigation  considered  a mean  flow  with  no  shear  in 
the  vertical.  Future  research  should  be  directed  at  investigating  the 
response  of  a vertically-sheared  flow  to  mountains.  Another  area  to 
be  examined  is  the  effect  of  condensation  and  the  release  of  latent  heat 
into  the  air  flow. 
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Figure  1 


ire 


a-  5 


level  2 


V/.T4- 
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The  sigma  (0)  coordinate  system  as  used  in.  the  two-layer 
version  of  the  model.  Also  shown  is  the  vertical  distri- 
bution of  the  variables. 
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Figure  10.  9°  x 1.5  km  mountain  range  case.  Six-layer  fine  resolution 

model  streamlines  for  a = .91667.  Dashed  lines  depict  out- 
line of  mountain  range  and  mountain  ridge  line. 
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Figure  11.  9°  x 1.5  km  mountain  range  case.  Six-layer  fine  resolution 

model  streamlines  for  a = .08333. 
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Figure  14.  9°  x 1.5  km  mountain  range  case.  Analytic  v profile  for 

a = 0.0  (dashed  curve)  vs.  six-layer  fine-resolution 
model  v profiles  for  a = .083333  (solid  curve)  and  a = . 
(dotted  curve) , all  at  44°N. 


Figure  17 „ 9°  x 1.5  km  mountain  range  case.  Six-layer  coarse 

resolution  model  streamlines  for  a = .75. 
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Figure  20.  3°  x 1.5  km  mountain  range  case.  Six-layer  fine  resolution 

model  streamlines  for  o = .91667. 


Figure  24.  3°  x 1.5  km  mountain  range  case.  Streamlines  at  0 

from:  a.  six-layer  fine  resolution  model 

b.  two- layer  fine  resolution  model 
(dashed  lines  indicate  the  topography) . 
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